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Abstract 

Incoherence in the controUed Hamiltonian is an important Hmitation on the precision of coherent control 
in quantum information processing. Incoherence can typically be modelled as a distribution of unitary 
processes arising from slowly varying experimental parameters. We show how it introduces artifacts in 
quantum process tomography and we explain how the resulting estimate of the superoperator may not be 
completely positive. We then go on to attack the inverse problem of extracting an effective distribution 
of unitaries that characterizes the incoherence via a perturbation theory analysis of the superoperator 
eigenvalue spectra. 

PACS numbers: 03.67.-a, 03.65.Yz 
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I. INTRODUCTION 



One of the biggest challenges in Quantum Information Processing (QIP) is the precise control 
of quantum systems. Errors in the control are conveniently classified as coherent, decoherent, and 
incoherent Q]. Coherent errors are systematic and differ from the desired operation by a unitary 
operation. Decoherent errors can be expressed by completely positive (CP) superoperators ^ and 



can be counteracted by techniques such as Quantum Error Correction (QEC) y, |^. An incoherent 
process can also be described by a completely positive superoperator [l|, |5|. The apparent non- 
unitary behavior of the incoherent process arises due to a distribution over external experimental 
parameters. The incoherent process is described by a superoperator S acting on Liouville space 
which can be written, when acting on columnized density matrices \p) obtained by stacking their 



columns on top of each other from left to right 



, as 



p{a)U{a) (8> U{a)da, 



(1) 



where p{a) is a probability density, i.e. the fraction of quantum systems within an ensemble that 
sees a given U{a) within an interval da, J p{a)da = 1, and U denotes the complex conjugate of 
U. This decomposition of a CP map into unitary Kraus operators is sometimes called a random 
unitary decomposition (RUD) 0|. A RUD exists for an incoherent process, but such a decom- 
position is sometimes possible even for a very general decoherent process [a] when there is no 
correspondence between the unitary operators in the decomposition and some actual distribution 
of associated experimental control parameters a. The distinction between the two therefore is prac- 
tical, and depends primarily on the correlation time of the variation of experimental parameters. 
If the latter quantity is longer than the inverse of the typical modulation frequency, the process 
falls into the class of incoherent noise ^, The point of making this distinction is that, whilst 
correcting for decoherent errors requires the full power of QEC, in practice incoherent noise effects 
are often reduced directly through the design of the time-dependence of the control fields. This is 
possible since the operators underlying the incoherence U (a) are assumed to be time-independent 
over the length of the expectation value measurement. Common approaches for instance in Nu 
clear Magnetic Resonance (NMR) include composite and adiabatic pulses ^ 
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18i | on composite pulses finds a great 



Furthermore, the work done by Tycko |Q| and Jones 
application in QIP since the schemes proposed are universal and therefore work regarless of the 
input state. 

In this paper, we demonstrate how incoherent errors introduce particular limitations to Quan- 
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FIG. 1: Quantum process tomography method. Starting from thermal equihbrium, a complete set of input 
states {pin} are prepared using a set {Sin} of control sequences. The input and output states are measured 
using a set of readout pulses {Uro} to rotate the density matrix into observable components {Mi„} and 
{Mout}- The process to be probed represented by a superoperator Sexp is then applied to these input states 
to obtain a corresponding set of output states {pout}- The measured map Sobs is then computed by right 
multiplying the matrix of output states by the inverted one of input states (see text for further details) . 
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due to the correlations they introduce with an 



environment in the QPT input states. Prior work has been devoted to the study of the implica- 
tions of such correlations in the system's reduced dynamics 22, S^,!^]. However, to our knowledge 



they have not explicitly been studied within the context of QPT to show that the method may 
output non-completely positive (NCP) maps. If the existence or origin of such noise is unknown, 
it is shown in 1] by means of explicit examples how one can eventually infer qualitative informa- 
tion, say its symmetry, about the probability distribution underlying the incoherent process from 
superoperator eigenvalue spectra. In this paper we tackle the inverse problem of extracting an 
effective probability distribution p{a) representing the incoherent superoperator, given a model 
for the source of the incoherent noise in the system. Such information is crucial for counteracting 
the incoherent errors which, due to their slow variation would otherwise persist during the 
experiment and quickly accumulate. 



II. INCOHERENT NOISE AND QUANTUM PROCESS TOMOGRAPHY 

QPT measures the experimental map associated with the implementation of a desired quantum 
operation by passing a complete set of input states through the gate and measuring the corre- 
sponding output states (see Fig. 1). 

This procedure is important for the experimental study of noise processes and for the design 
of quantum error correcting codes P|. If incoherent noise is present in the preliminary step of 



3 



QPT, then the prepared (input) states will be classically correlated with the control parameter 
a characterizing the incoherence (Eq 1). Furthermore, if the correlation time of the noise in 
the control parameter is long compared to the coarse-grained time at which the evolution of the 



system is monitored, then the subsequent dynamics is non-Mar kovian j25|. More specifically, the 
gate applied to the input state (i.e., the gate being characterized by QPT) will be correlated with 
the same (slowly- varying) control parameter, and therefore also correlated with the input state to 
which it is applied. In such cases, the measured dynamics are not guaranteed to be completely 
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positive, and need not correspond to a linear map 

More generally, any correlations arising from non-Mar kovian dynamics, whether quantum or 
classical, can lead to incorrect interpretations of the data obtained from QPT. The "environment" 
which we assume from the outset is correlated with the system is in general defined by the degrees 
of freedom that are not part of the system. For example, the different spatial locations of individual 
qubits in a NMR ensemble, or the bosonic bath producing the fluctuations of the gate charge in a 
superconducting qubit [2^. 

The basic issues can be seen by exploring the QPT of a spin-1/2 system coupled to a second 
spin- 1/2 system B as its environment. Borrowing the example used in the initial density 
matrix of the total system may be written as 

PAB = ^i^AB + aiUi 0Ib + l3jlA ® CFj + lijCTi ® aj), (2) 

where / and fij denote the identity and Pauli matrices respectively. The density matrix of the 
system A is then obtained by tracing over the environment B: 

PA = ^ilA + o:iai). (3) 

The dynamics of the whole system pab is a unitary evolution Uab, so that the resulting density 
matrix of A is 

p'a = ^{p\Uab{pa® Pb)U\i^\p) 

where 7^ = (7^^ - aif5j)/A and ps = Tta{pab)- Let ps = YjvPvW){i'\ and M^y = {p\Uab\v). 
Then the above expression can be reexpressed as 



+ ^{l^\UABlijCFi (g) crjU\j^\n) 

The first line therefore corresponds to the Kraus operator sum form [2^ of the evolution when 
initial correlations are not present, while the second line represents the contribution from these 
correlations. We can easily see that there exist Uab, e.g., the swap gate, for which these initial 
correlations are not observable. In general, however, initial correlations cause the map to be non- 
linear or NCP. 

Within the context of QPT, we now investigate an explicit example of how NCP superoperators 
can arise. We take a set of 4 initial density matrices pab such that /O^ is the same in each case, 
and where the input states pA span the Hilbert space of the system A (as required by the QPT 
procedure) 

Pab = {lAB + PI®a,)/^^ p\i^ = I/2 

Pab = {Iab + aa^ ® I + 131 ® + -fa^ ® (yz)l^ =^ Pam = + aa^:)/! 
p\b = {Iab + aay (g) I + (31 (g) + -fay (g) a^)/4: =^ p\ i,^ = {I + aay)/2 
Pab = i^AB + aa^ I + (31 ^ + 'jdz ® (Tz)!^ =^ PA,in = + ao-^)/2 



where in each case pB = {I + (iaz)/2. With the example Uab = 6^*4'^^®°"^, the corresponding 4 
outputs are p\ = 1/2, p\ = {I + ^(Ty)/2, p\ = {I — ^ax)/2 and p\ = {I + aaz)/2. We write 
the density matrices p as vectors in Liouville space in the Zeeman basis j^, which are obtained 
by first writing the density matrix in the Zeeman basis and then stacking their columns on top of 
each other from left to right. We refer to the resulting vector simply as the "columnized density 
matrix", and will denote it as a ket \p). If we set a = /3 = 0.5 and 7 = 0.6, the map S is 



Output 1 Output 2 Output 3 Output 4 

0.5 0.5 0.5 0.75 



0.3i 
-0.3i 
0.5 0.5 



-0.3 
-0.3 
0.5 






0.25 



S- 



Input 1 Input 2 Input 3 Input 4 

- 0.5 0.5 0.5 0.75 

0.25 0.25* 

0.25 -0.25i 

L 0.5 0.5 0.5 0.25 



S 



10 

l.2i 

-1.2j 

1 



(4) 



which is in general non-linear, since one can no longer predict the output for an arbitrary input 
state given the action of the gate on these four specific input states alone. However, the map can 
be considered to act linearly on the system A Hilbert space, i.e. on the linear combinations of input 
states which contain the right correlations with the environment. For instance, the action of the 
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gate on the input state (/9^+p^)/2 = {I+^{ax+(Ty)) /2 can be computed by using the above matrix 
expression if the total input state is pab = {Iab + ^i^^x + (Jy) ® I + {31 ® az + + ay) ® Uz) / ^■ 
This result conveniently allows one to treat the map as linear. If treated as linear, the Choi matrix 
[3| C = J^i'j^oi-^ij ^ ^)'^{^ -^ij)) where is the dimension of the system's Hilbert space and 
Eij is the N X N elementary matrix (with a "1" in the ij — th position and zeros elsewhere), 
corresponding to the superoperator S is not positive semidefinite and consequently S can not be 
CP U. 

It is suggested in [y| how the NCP part of the superoperator can be removed, namely, by 
removing the negative eigenvalues of the Choi matrix and then renormalizing so that the trace 
is equal to the dimension of the Hilbert space, N. We shall call this method CP-filtering. The 
Choi matrix corresponding to S* in Q has two non-zero eigenvalues (2.2,-0.2). In this exam- 
ple, the CP-filtering procedure replaces the negative eigenvalue by and then renormalizes the 
new Choi matrix to trace N = 2. The CP-filtering method outputs one unitary Kraus operator. 
On the other hand, if no initial correlations were present, the superoperator would be equal to 
S = diag(l, 0.5i, — 0.5i, 1), and the corresponding Choi matrix would have two positive eigenval- 
ues (0.5,1.5), yielding two Kraus operators. This superoperator is therefore completely positive. 
Unsurprisingly, we thus observe that the superoperator obtained from the CP-filtering procedure 
is not equal to the superoperator obtained by removing the initial correlations. Therefore, unless 
the initial correlations are very small, CP-filtering is a fairly uncontrolled procedure, giving CP 
superoperators that may significantly misrepresent the true quantum dynamics. 

We now take the case a = /3 = 7 = 0.5 with the same Uab as previously. Although initial 
correlations are still present, the superoperator obtained by the above QPT method, without CP- 
filtering, is CP with one Kraus operator. In contrast, if the correlations in the initial states are 
removed whilst keeping all other things equal, the superoperator obtained is CP with two Kraus 
operators, not one. So even when initial correlations are present, the superoperator obtained via 
QPT may be completely positive. Therefore, one cannot rule out the presence of initial correlations 
merely by the existence of a valid Kraus operator sum form via QPT. Initial correlations can 
masquerade as CP maps, and in reality the process may not be linear with respect to arbitrary 
input states. 

To summarize, the results of this section are: (i) incoherent errors introduce correlations between 
the system and the environment in the QPT input states which can persist during the implemented 
transformation, (ii) these correlations can yield non-completely positive superoperators or non- 
linear maps, (iii) the CP-filtering method suggested in |a| is not equivalent to removing these 
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TABLE I: Summary of the two examples given in section |n] CPF is abbreviation for CP-Filtering which is 
an algorithm for converting a non CP superoperator to a CP one. "Corr" is abbreviation for correlations 
(between subsystems A and B). 

initial correlations, and (iiii) initial correlations can masquerade as CP maps which misrepresent 
what is in reality a non- linear process with respect to the input states. See Tabled for a summary. 
This simple analysis explains the apparent NCP behavior measured in experiments reported in 
20, 



I29I 1. This motivates the need to characterize the incoherent noise, and to provide ways to 
correctly interpret QPT data. If the noise can be successfully characterized, we may use this 
information to better counteract the noise in the first place. 



III. EXTRACTING PROBABILITY DISTRIBUTIONS FROM SUPEROPERATOR 
EIGENVALUE SPECTRA 



By applying a first order perturbation theory analysis of the eigenvalues of superoperators, 
we now present a method to extract the probability distribution profile p(a) of unitary matrices 
present in incoherent processes. For an incoherent noise to be refocused [lO|, knowledge about the 
noise is a priori required. If qualitative information about the inhomogeneity in the Hamiltonian 
is known, spectroscopic techniques can be used to obtain the missing quantitative information. In 
the following analysis, we assume that the physical origin of the incoherent noise is unknown or 
hidden due to the complexity of the system-apparatus interactions, but that a mathematical model 
is presumed. 

As presented in the introduction, an incoherent process implies a random unitary distribution. 
Incoherent processes are thus unital, which means that the maximally mixed density matrix is 
left unchanged. A linear, completely positive, trace preserving and unital map is called a doubly 
stochastic map [s^. Although a single necessary and sufficient condition for a doubly stochastic 
map to possess a RUD has not been found to our knowledge, examples of doubly stochastic maps 
which do not possess a RUD are reported in PiQl- However, since many decoherent unital processes 
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can be modelled by a stochastic Hamiltonian, i.e. semiclassically, we believe that many instances of 
decoherent processes can have a RUD. This belief is supported by the following two facts. Any two 
density operators {p,p') connected by a doubly stochastic map, p' = A(p), can always be related 
by a transformation of the form p' = YliPi^iP^i ■< where YliPi — 1 ^-^^ {Ui} is a set of unitary 
operators js^. Furthermore, all unital maps for a two-level quantum system always have a RUD 



A. Perturbation Theory Analysis of the Eigenvalue Spectra 



In what follows, we take examples from NMR 



3ll | physics where the main source of incoherence 



comes from Radio Frequency (RF) power inhomogeneity. Due to the s|)atial extent of the sample, 
individual spins during the course of a RF field see different powers 1] and evolve according to 
different unitary evolutions with different characteristic frequencies. Note that identical spins can 
have different resonance frequencies due to inhomogeneity in the static magnetic field within the 
ensemble, which is another source of incoherence. However, as shown in _^], the non-unitary 
features arising from this static external field inhomogeneity are usually much smaller than those 
arising from RF inhomogeneity and will be therefore ignored in this example. Finally while the 
distribution of RF fields can be easily measured via a nutation experiment on a single spin, the 
method presented here is quite general and can for instance account for the correlation between 
multiple sources of incoherence (several RF fields, DC field etc.). 

Let n be the number of spin-1/2 particles in the ensemble, and Uk denote the unitary operator 
generated by the RF field in the kih. frequency interval of the RF amplitude profile. The eigenvalues 
of the superoperator S = Uk ® are products of the eigenvalues of with those oi Uk- This 
yields 2" eigenvalues that are equal to unity and (2^"~^ — 2"~^) pairs of eigenvalues (A, A). In 
general, the eigenvalues of CP superoperators come in conjugate pairs, but only in the case of 
unitary superoperators do all the eigenvalues lie on the complex unit circle. 

The incoherent process resulting from an inhomogeneous distribution Uk processes is given 
by the superoperator S = J2kPk^k ® Uk, where pk is the fraction of spins that sees the unitary 
evolution Uk- The more broadly the {pk} are distributed, the larger the degree of inhomogeneity 
in the evolution, and the more incoherent noise enters into the evolution. Estimates of the actual 
eigenvalues of = Y2k Pk^ k®Uk will now be obtained using non-degenerate first-order perturbation 
theory. Because the RF pulses are not perfect, even in the absence of RF field inhomogeneity, 
we may assume that the unperturbed eigenvalues are generically non-degenerate. The unitary 
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operator Uk generated by the RF field acting at position k may be written in exponential form as 
Uk = e~*^* * where represents the effective Hamiltonian of the evolution over the period t for 
which the pulse is applied {H has been set equal to 1). Defining Hq to be the unperturbed (and 
desired) Hamiltonian, the eigenvalues and eigenstates \4>j) of Hq satisfy the eigenvalue equation 



Uo\^j)=e-"f'^%), (5) 

where Uq = ex.p{—iHQt). The Hamiltonian of a particular Uk is assumed to be a perturbation of 
the desired, homogeneous Hamiltonian 

Hk = Ho + Kk, (6) 
where is the perturbation. To first order, the new eigenvalues of Hk are 

^j,fe = + (7) 
and the corresponding eigenvalues oi Uk are 

Uk\hk) = e-'^'''%,k). (8) 
To first order, the spectral decomposition of S is 

S=J2Pk (e'^'"'' '\<Pm){<Pm\ ® 6-^^-'=* 10,-) (</),■ l) (9) 
k,m,j 

and the eigenvalues of S are given approximately by 

k 

k 

We now imagine the scenario where Kk is given by K^t = K. This result would in fact be 

exact for one spin on resonance. In this case, is the parameter a defined in the introduction 

(which parameterizes the inhomogeneity) and represents the normalized RF power deviation from 
the desired power loq. Defining Ao; = and Kj^ = {<^j\K\(pj) — {(t>m\K\(t>m) , the previous 

equation in the continuous limit becomes 

= e-i{4>j-4>m)t J p(A(j)e-'^^-^^dAu;. (10) 

We see in this case that to first order the eigenvalue Xjm is just the unperturbed eigenvalue 
e-i{'t>j-'t>m)t times the Fourier transform of the RF distribution profile evaluated at Kjj^. This 



result demonstrates that the probabihty distribution profile p(Aw) of an incoherent process can 
be determined, within some degree of approximation, from the eigenvalue structure {Xjm} of an 
experimental superoperator, given some model for the incoherence K. Knowing K would indeed 
allow one to build the correspondence between Xjm and Kjm, and then to determine p{Alo) by 
performing an inverse Fourier transform. Of course, this result holds for general K only when the 
perturbation is in the first order regime, and when the unperturbed eigenvalues are non-degenerate. 
But if K approximately commutes with Hq, then the first-order perturbation in the eigenvalues 
is close to an exact correction, and the above analysis gives a very accurate description of the 
incoherent process. 



B. Recovery of the Profile 

We now demonstrate via a numerical example how one can recover the profile p{Auj) from the 
eigenvalue spectrum of a measured superoperator. In the theory derived in the previous subsection, 
Xjm is the data, i.e., the eigenvalues from the measured superoperator. A model K is needed for 
the perturbation, while (j)j and \(j)j) are known through the knowledge of Hq. Formally solving for 
p{Alo) from ((TUIl . 

p{Auj) = Xjme'^^^-^^-^'e'^'^-^^^dKjm. (H) 

The different eigenvalues Xjm multiplied by g^^'^^"'^'"-'* therefore allow us to obtain the complex 
function of Kjm which we shall call /, corresponding to the Fourier transform of the distribution 
profile. 

For the numerical demonstration of this technique we take a 3-qubit system. We choose Hq and 
K such that \{(i)]\K\(l)„) / — (j)n)\ — 0.1 for n ^ I and [Hq, ~ so that first order perturbation 



theory can be used 



32l |. We then use a measured RF inhomogeneity profile shown in Fig. 2 to 



construct the following superoperator acting on Liouville space: 

S = Y^ p{Auj)U{Auj) U{Alo) . (12) 

Provided with this superoperator, we compute its eigenvalues and plot them on the Argand 
diagram (see Fig. 3). The perturbation K is such that the first order limit condition is fulfilled, 
but that the size of its diagonal elements in the Hq eigenvectors basis is large enough to generate 
significant dephasing and attenuation in the eigenvalues. The correspondence between Xjm and 
Kjm is needed to recover the distribution profile (fTT|) . The correspondence can be established by 
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FIG. 2: Radio Frequency inhomogeneity profile used to construct the superoperator S {J p{ALo)dALU = 1). 
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FIG. 3: Eigenvalue spectra of the incoherent process 5* and the desired Sq. The dots are the eigenvalues 
of the incoherent superoperator S and the diamonds the ones of the desired unitary superoperator 5*0 (on 
the unit circle). Also shown is the correspondence between the miperturbed and perturbed eigenvalues. 
The perturbation K is small enough so that the first order limit condition is fulfilled but large enough to 
substantially dephase and attenuate the eigenvalues. 

first computing {(pj, (j)m\S\(j)j, (pm), then searching for the eigenvalue of S closest to it. This allows 
us to make the correspondence between one unperturbed eigenvalue with eigenvector \(j)j,(pm) = 
\4>j) ^ \4>m) (obtained from the knowledge of Hq) and one eigenvalue of S. The function / = 
J p{ALo)e~^^'^^^"'dAuj with respect to Kj^n can then be constructed. The real and imaginary parts 
of that function are plotted in Fig. 4. Note that we ignore the degenerate points at / = 1 because 
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FIG. 4: Plot of / = / p{ALu)e 



with respect to Kj„i (real and imaginary parts). The function is 



conjugate symmetric with respect to as expected, so that its inverse Fourier transform, which should be 
a probability distribution, is real. The point at Kjm = was added to avoid a DC offset in the reciprocal 
Fourier domain. 

they do not provide any information about p{Auj) other than normalization. The 64 eigenvalues 
of the superoperator minus the 8 degenerate ones equal to 1 (at Kjm = 0), plus 1 eigenvalue 
added at Kjm = to avoid a DC offset in the reciprocal domain, yield a complex function / 
made of 57 unequally spaced points. The function / is conjugate symmetric with respect to 0, 
which is consistent with the fact that we are supposed to recover a probability distribution, i.e. 
a real function, after computing the inverse Fourier transform. To perform the inverse Fourier 
transform of a function sampled at unequally spaced points, we used an algorithm prescribed in 
[3^ . The result is shown in Fig .5. The width of the probability distribution and its skewness are 
recovered to a good extent, the discrepancy being due to the lack of information about the function 
/. It is worth mentioning that with 57 sample points, the window of Kjm values should be large 
enough to allow low frequency components of the profile to be reliably extracted. If the incoherent 
perturbations were very small, there will not be as many large values of Kjm, and therefore less 
low frequency information would be available. However, the perturbations can be made larger 
without changing the mathematical model, by simply repeating the control sequence several times, 
provided other noise mechanisms do not play a significant role. In addition, more points could be 
used to get a better sampling resolution, with, for instance, a 4-qubit superoperator yielding 241 
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FIG. 5: Inverse Fourier transform of the function plotted in Fig. 4 shown together with points from Fig. 2. 
The width of the profile in addition to its skewness are recovered to a good extent. 

points. One cannot have arbitrarily many points, however, because the correspondence between 
the eigenvalues of S and those of Sq could quickly become impossible to establish unless a very 
good knowledge of the perturbation K is available. The density of points in the Argand diagram 
becomes so large that eigenvalues can easily become confused. Here, the 3-qubit superoperator 
was enough to recover the essential features of the probability distribution. 

If Hq is not exactly known, and in fact a constant offset Hamiltonian which is proportional to 
K is present, then a different function / is obtained : 



where /3 is a constant real number. Taking the inverse Fourier transform of / would reveal a 
distribution p{Auj) centered around /3 rather than 0, indicating that Ho + f3K is in fact the unper- 
turbed Hamiltonian. Perfect knowledge about Hq as a result is not required provided the offset is 
approximately proportional to K. 

It is important for this method to work that the model K chosen a priori is a reasonably 
faithful one, and that it approximately commutes with Hq. This ensures that the applied first- 
order perturbation theory is valid. In the extreme case, in which K anticommutes with Hq, 
{(l)j\K\(f>j) = and no "data" would be available for analysis. It is also worth mentioning 
that this method is, needless to say, not scalable. However, in many instances, as in NMR, the 
inhomogeneity features are apparatus dependent, so that reasonably small physical systems can 
be used to probe them. The scalability of the method is not necessarily a requirement. The 
distribution of some control parameters, once obtained, can be valuable in designing robust control 




(13) 



13 



sequences [jj] for larger and more complex systems. 



IV. CONCLUSION 

Here we reviewed that when incoherence is present during the preparation of the input states for 
QPT, the resulting correlations between the system and the "environment" can play an important 
role on the subsequent system's dynamics. The map obtained by right multiplying the matrix 
of output states by the inversion of the matrix of input states still has a meaning, but a correct 
interpretation of the measured data (or transformation) requires an analysis of the incoherence 
effects affecting the tomographic procedure. In particular, the measured map needs not be CP. If 
quantitative information is missing, our perturbation theory analysis of superoperator eigenvalue 
spectra can be used to determine an effective distribution of unitaries characterizing the process, 
provided a good mathematical model is available. While this requires a significant effort to measure 
a 3-qubit superoperator, it is certainly feasible within present experimental capabilities. Lastly, the 
knowledge of the distribution of control parameters should finally allow us to design more efficient 
control sequences aimed at counteracting these deleterious effects. 
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